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We have carried out an investigation of the environments of low redshift HII 
galaxies by cross-correlating their positions on the sky with those of faint field galaxies 
in the Automatic Plate Measuring Machine catalogues. We address the question of 
whether violent star formation in HII galaxies is induced by low mass companions 
by statistically estimating the mean space density of galaxies around them. We argue 
that even if low mass companions were mainly intergalactic HI clouds, their optical 
counterparts should be detectable at faint limits of the APM scans. 

A significantly positive signal is detected for the HII galaxy- APM galaxy angular 
cross-correlation function, but the amplitude is poorly determined. The projected 
cross-correlation function has higher signal-to-noise, and suggests that the amplitude 
is slightly lower than for normal field galaxies. This implies that these bursting dwarf 
galaxies inhabit slightly lower density environments to those of normal faint field 
galaxies, consistent with other studies of emission line galaxies. This suggests that 
in these dwarf starburst galaxies, star formation is not triggered by unusually strong 
tidal interactions, and may have a different origin. 

Key words: galaxies: statistics - galaxies: starburst - galaxies: clustering - galaxies: 
environment 
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1 INTRODUCTION 



The starburst phenomenon is observed in a large num- 
ber of extragalactic objects from giant HII regions in irregu- 
lar galaxies and late type spirals, to galaxies entirely domi- 
nated by the violent massive star formation region as in the 
case of Starburst Galaxies. The class of starburst galaxies 
comprises a large range in luminosity, mass, heavy element 
and dust content, as well as morphology. Classical starburst 
or nuclear starburst galaxies typically have an intense re- 
gion of violent star formation in the nucleus of an otherwise 
normal spiral galaxy (Balzano 1983). At the high luminos- 
ity end of starburst galaxies, ultraluminous IRAS galaxies 
(Soifer et al. 1987) are strongly interacting giant systems 
(Melnick & Mirabel 1990) where most of the radiation is 
emitted in the far-infrared due to reprocessed UV radiation 
by the large content of dust particles. HII galaxies, on the 
other hand, are dwarf galaxies in a bursting phase of star 
formation of low luminosity and mass, low heavy element 
abundance and low dust content where the triggering mech- 



anism of the present episode of violent star formation is not 
as obvious (Telles & Terlevich 1995). 

Earlier searches for bright galaxies near to HII galaxies 
found a deficit of L > L* galaxies within 1 Mpc (Campos- 
Aguilar & Moles 1991, Campos- Aguilar et al. 1993, Vilchez 
1995, Pustilhik et al. 1995; Telles & Terlevich,1995). HII 
galaxies are not associated with giant galaxies, therefore 
they are not tidal debris of strongly interacting systems. 
Auto-correlation analyses of strong emission line galaxies 
(Iovino, Melnick & Shaver 1988, Rosenberg et al. 1994, Love- 
day et al. 1999) find a low clustering amplitude, suggesting 
that HII galaxies tend to populate regions of low galactic 
density. 

Their optical properties are dominated by the massive 
star forming region, as shown by their strong emission line 
spectra superposed on a weak blue continuum. The proper- 
ties of the underlying galaxies in these systems are similar to 
late type dwarf galaxies such as dwarf irregulars or low sur- 
face brightness dwarfs (Telles & Terlevich 1997). The most 
luminous HII galaxies, classified as Type Ps by Telles, Mel- 
nick and Terlevich (1997), show signs of disturbed morphol- 
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Figure 1. The distribution of HII galaxies and APM data on the sky in an equal area projection in equatorial coordinates. The solid 
circles are the positions of the HII galaxies. The square fields are the APM scans of the UK schmidt plates. The dashed lines show lines 
of constant RA and DEC. 



ogy such as distorted outer isophotes, tails or irregular fuzz, 
while the low luminosity Type IPs are regular and compact. 
Although there is a clear case for a morphology-luminosity 
relation, neither type of HII galaxies shows conspicuous ev- 
idence of bright companions in their neighbourhood. The 
few HII galaxies found to have a bright neighbour (maybe 
by chance) are all of Type IPs of regular morphology, con- 
trary to what one would expect if interactions caused the 
morphological disturbances as seen in Type I's (Telles & 
Terlevich 1995). 

A popular hypothesis is that interactions between 
galaxies are the triggers of starbursts and they may also 
cause the current burst of star formation in HII galaxies. 
Giant galaxies, however, are not found in the immediate 
vicinity of HII galaxies, thus are improbable candidates for 
triggering agents. A possible and appealing alternative was 
presented by Melnick (1987) who proposed that high res- 
olution 21cm maps were needed to investigate the role of 
collisions between intergalactic neutral hydrogen clouds in 
the formation of these objects. Brinks (1990) also hypothe- 
sized that other dwarfs or intergalactic HI clouds could be 
the triggering agents. Taylor et al. (1995, 1996) using the 
VLA detected 12 HI companions around 21 HII galaxies, 
while only 4 Hi-rich companions were detected around a 
control sample of 17 quiescent low surface brightness dwarfs 
(Taylor 1997). As also pointed out by these authors, some 
questions remain intriguing from this: Why are these 9 out 
of the 21 HII galaxies with no companions violent forming 
stars now ('bursting')? Why are these 4 out of the 17 LSBGs 
with companions not 'bursting' ? 

HI surveys find that all the HI detections have an optical 
counterpart. That is, all the sources found in 21 cm surveys 
are nothing else than normal galaxies (c.f. Zwaan et al. 1997, 



Zwaan, 1999, Hoffman, 1999). No free floating intergalactic 
HI clouds were detected in such surveys. Thus, we have car- 
ried out a further investigation of the galaxy environments 
of a unbiased sample of over 160 low redshift HII galax- 
ies by cross-correlating their accurate position in the sky to 
faint field galaxies (15 < bj < 20.5) in the Automatic Plate 
Measuring Machine (APM) galaxy catalogue. For the mean 
redshift of our HII galaxy sample we detect galaxies down 
to Mb ~ —14.5. Using the relation between optical magni- 
tude Ms and HI mass Mm for late type galaxies, given by 
Rao & Briggs (1993) [log X hi M q = 2.72 - 0.36M B ], we es- 
timate that we are not missing any cloud with mass greater 
than 10 s Mq. This is comparable to the lower limits of HI 
companions found by Taylor et al. (1995, 1996). Hence our 
present study should detect any possible low mass candi- 
dates to act as tidal triggers. 

In Section [] we describe in more detail the data-sets 
used in the present work and we present the details of our 
calculations in Section |3| Finally, we show some of our con- 
clusions in Section m. 



2 DATA 

2.1 HII Galaxy Sample 

The HII galaxy sample used in this paper is taken from 
the Spectrophotometric Catalogue of HII Galaxies (SCHG; 
Terlevich et al. 1991). Most of the objects in this database 
have been selected from the Tololo survey (Smith, Aguirre 
& Zemelman 1976), and the University of Michigan survey 
(MacAlpine & Williams 1981). The catalogue also contains 
a number of objects which are not classified as Seyfert galax- 
ies selected from the Markarian list of galaxies with strong 
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2.2 APM galaxy sample 




Figure 2. The distribution of redshifts for the HII galaxies (his- 
togram) and the model distribution (equation |l] used to calculate 
the model correlation function (dashed line). 



ultraviolet continuum (Markarian, Lipovetskii & Stepanyan 
1981 and references therein), as well as some blue objects of 
Zwicky's catalogue of compact galaxies (Zwicky 1971). The 
total catalogue contains spectra of over 400 emission line ob- 
jects found in objective prism surveys using Illa-J emulsion 
through their [OIII]AA4959,5007 and/or [OII]AA3726,3729 
lines. From these, about 300 are HII galaxies. The remain- 
der are Giant HII regions, and Starburst nuclei or emission 
line objects classified as Seyfert nuclei from their position in 
the emission-line-ratio diagnostic diagrams. 

Most of the objects in this sample cover only two specific 
regions of the sky. For instance, the Michigan survey covers 
a 10° band around the celestial equator, while the Tololo 
survey concentrates in the region —27° < S < —43°. This 
is illustrated in Figure [j] which shows the distribution of 
our HII galaxies in the sky. For the present study we ended 
up with 163 HII galaxy centres for which there are APM 
scanned UK Schmidt plates. The actual centres for the HII 
galaxy were carefully identified for each APM field, thus 
assuring that the HII galaxy is not counted as a companion 
of itself. 

The redshift distribution of these galaxies is plotted in 
Figure ^, which also shows the best-fit model redshift dis- 
tribution of the form 



N(z) oc — exp 



z_\ 



(1) 



The mean redshift is 0.03. The typical absolute magnitude 
is Mb ~ —18. Throughout this paper we use the current 
value of H = 65 kms" 1 Mpc _1 (Suntzeff et al. 1998). 



Our sample of faint field galaxies was selected from the APM 
Galaxy Survey, which is described in detail by Maddox et 
al. (1990). The sky covered by the APM galaxy survey has 
been extended since the original publication of the survey: 
the south galactic pole part of the survey now covers a solid 
angle of 6250 square degrees, and is based on 269 UKS J 
plates scanned by the APM machine; the north galactic cap 
area covers 750 square degrees from scans of 30 UKS plates 
centred with 9 h < a < 15 h and -5° < 6 < 0°. The fields 
covered are shown by the dotted squares on Figure 

The data consist of positions accurate to < 1" and mag- 
nitudes accurate to ~ 0.1 mag for over 50 million images 
brighter than a magnitude limit of bj = 21.9. The galaxy 
photometry has been corrected for several systematic ef- 
fects and have no detectable systematic errors more than 
~ 0.04 magnitudes rms. The galaxy sample selected from 
the survey data at a magnitude limit of bj — 20.5 has a 
completeness ~ 90-95%, stellar contamination ~ 5% and 
negligible dependence of the galaxy surface density on dec- 
lination or galactic latitude (Maddox et al 1996). At this 
magnitude limit the redshift distribution is well described 
by equation |l| with the mean z = 0.16 (Maddox et al 1996). 

We extracted APM measurements for a 2° square area 
around each of the HII galaxies. For the central 10' square we 
visually cross-checked the APM catalogue list against images 
from the Digitized Sky Survey (DSS), and rejected multiple 
detections and noise images. This provides a reliable galaxy 
catalogue around each HII galaxy. 



3 ANALYSIS 

3.1 The angular cross-correlation function 

We have measured Wh g {9), the angular cross-correlation 
function between the HII galaxies and an apparent magni- 
tude limited sample of neighbouring galaxies. We estimated 
Wh g by counting the number of galaxies Nh g as a function 
of angular radius 9 from the central HII galaxy, and com- 
paring this to the number Nh r counted for a catalogue of 
uniform random postions. We used ten times as many ran- 
dom points as galaxies in order to reduce their contribution 
to the counting errors, and then rescaled the count down by 
a factor ten. The cross-correlation function is then given by 



W hg(0) 



N hr (6) 



- 1 



(2) 



We also used the simpler direct estimate using the mean 
surface density of field galaxies, A/ , and the area of each 
annulus, 



N hg (e) 



n cm Nit{{0 - A) 2 



'-)) 



(3) 



where n cen is the number of HII galaxies used as centres. 

The main results are shown in Figure []. The filled points 
show w from equation |^ and the open points from the direct 
estimator, equation^. This gave essentially indistinguishable 
results, showing that there are no significant systematic bi- 
ases in our sample. The error bars are estimated from the 
Poisson noise in each radial bin, scaled by the integral of w, 
e w = (1 + 2-RNJ 2 {9))/^/N h ~ g , where J 2 = /„" 9w{9)d9. This 
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Figure 3. The angular cross-correlation between the HII galaxies and the faint APM field galaxies. The filled points show the results 
using equation ^, and the open points are from equation |i| The error bars are estimated from e w = (1 + 2nAfJ2(9))/ N^ g . The line is 
the cross-correlation function predicted assuming that the HII galaxies are a random subsample of APM galaxies. 



is analogous to the error estimate for the spatial correla- 
tion function £ suggested by Kaiser (1986), and is a reason- 
able approximation for a weakly clustered distribution (see 
Hamilton 1993 for an extensive analysis of errors in corre- 
lation functions). It can be seen that Whg is significantly 
positive for angles 9 < 30', corresponding to an excess of 
galaxies near the HII galaxy positions compared to a uni- 
form distribution. This would be expected for any sample of 
galaxies, since we know that galaxies are clustered. 

We have calculated the expected cross-correlation func- 
tion assuming that HII galaxies cluster in the same way as 
normal galaxies, and this is shown by the line in Figure || 
This prediction is based on the APM correlation function 
for bj = 20, which is well fit by a power law at small scales, 
w(9) = A9 1 ~~ f , with 7 = 1.699 and Aapm = 0.0284 (Mad- 
dox et al 1996). We scaled the amplitude by a factor calcu- 
lated numerically from Limbers equation (Peebles 1980) 

(os _ /o°° IT N ^) N Az 3 )qr hg )(£) h (£) g dz h dz g 
WhA) [J °°N h ( Z )dz][J °°N g (z)dz] 

where x is the comoving coordinate at redshift z, Nh(zh) is 
the redshift distribution of the HII galaxy sample, N g (z g ) 
is the redshift distribution of the APM galaxy sample and 
£i r hg) is the spatial correlation function at separation r^ 3 . 
The two redshift distributions are given by equation ^ with 
z c chosen to match the observed distributions for HII galax- 
ies and the APM galaxies at b.j = 20.5. This gave the 
predicted amplitude Ah g ~ 0.01399. We also included a 



constant offset to correct for the integral constraint within 
the 2° fields (Smith, Boyle &Maddox 1995), giving a final 
w(6) = 0.01399(6»-°' 699 - 1.5), where 9 is in degrees. 

Figure [] shows that the estimated Whg (9) for HII galax- 
ies is consistent with the assumption that HII galaxies are 
clustered in the same way as normal field galaxies. The un- 
certainty in the estimated amplitude is rather large, but we 
can rule out the suggestion that HII galaxies are anticorre- 
lated with field galaxies. It is also clear that there is not a 
large excess of near neighbours around the HII galaxies com- 
pared to normal galaxies. This appears incompatible with 
suggestion that HII galaxies are triggered by tidal interac- 
tions with nearby low-mass galaxies. 



3.2 The projected cross-correlation function 

Since we know the redshift to each HII galaxy, we can es- 
timate the correlation function using the projected physical 
separation, a. The resulting projected cross-correlation func- 
tion, , is an integral over the spatial correlation function 

/oo 
i hg ({x 2 + a 2 ) 1 ' 2 ) dx (5) 
- oo 

We estimate Hh s using a method similar to that described 
by Saunders et al (1992). For each HII galaxy we count the 
excess neighbours compared to a random distribution using 
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Figure 4. The projected cross-correlation between the HII galaxies and the faint APM field galaxies, B,f lg (a)/a. The points show our 
measurements, with error bars estimated from the scatter between centres. The solid, dotted and dashed lines show the predicted 3 for 
ro = 5.1,3.7 and 2.7h~ 1 Mpc respectively. 



This is related to H 



E(a) 



J n(x)dx 
n(y) 



X(a) 



(6) 



(T) 



where n(x) is the average number of field galaxies per stera- 
dian per Mpc along the line of sight in the background cat- 
alogue at distance x, and the distance to the central HII 
galaxy is y. The approximation is essentially the small-angle 
approximation, but also involves several subtleties, as dis- 
cussed by Saunders et al (1992). Note that the different dis- 
tance to each HII galaxy means that relation between a and 
9 is different for each centre, and also the l/n(y) leads to a 
different weighting of the pair count from each centre. This 
means that S,h g is not simply a rescaling of Wh g - 

If £h g is a power law in r, £h g = (r/ro) -1 , then the pro- 
jected correlation function is given by H>, 9 (cr) = (ct/cto) 1 " 7 
where a 7 " 1 = r 7 r(i)r(^i)/r(^)! Hence 



E hg (a)/a 



:& 9 (a)r(i)r(2=i)/r(f). 



(8) 



Our measurement of ~Eh g {<j)/a is shown by the points 
in Figure ^. It is positive for cr < IMpc, showing that HII 
galaxies have more neighbouring galaxies than a uniform 
distribution: they are not isolated systems. As discussed in 
Section 3.1, we expect any sample of galaxies to have more 
neighbours than a random distribution, so we have calcu- 
lated the expected S/cr assuming that HII galaxies are clus- 



tered in the same way as normal galaxies. Our prediction is 
given by equation ^| with Chs( r ) = OV^l/i - 1 ) -1 " 7 ' L , which 
is a good approximation to the APM correlation function 
on small scales (Maddox et al. 1996). We have also sub- 
tracted a constant from the power-law to allow for the in- 
tegral constraint which applies to the data points. This is 
shown as the solid line in Figure [|. The observed values 
are consistent, although slightly lower than the predictions. 
This slight discrepancy can be interpreted as a reflection of 
the fact that the amplitude of the autocorrelation function 
for HII galaxies is lower than for normal galaxies. Iovino 
Melnik and Shaver (1988) find r = 2.7/i _1 Mpc for the HII 
galaxy sample used here; the corresponding E;,h is shown 
as the dashed line in Figure ^ If HII galaxies and normal 
galaxies follow the same underlying mass distribution with 
differing bias levels, the cross-correlation function should be 
the geometric mean of the two autocorrelation functions, 
so ^ hg = y/^hhigg- The dotted line in Figure [j| shows the 
equivalent Sh g - 

The lower amplitude emphasises the point made in Sec- 
tion 3.1, that there is not a large excess of near neighbours 
around the HII galaxies compared to normal galaxies. Again 
this is incompatible with suggestion that HII galaxies are 
triggered by tidal interactions with nearby low-mass galax- 
ies. 
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4 CONCLUSIONS 

Our main results are: 

(i) Both the angular and projected correlation functions 
are significantly positive, so HII galaxies are significantly 
clustered. This is what you expect to find for any sample of 
galaxies. 

(ii) The angular measurements have large uncertainties, 
but are consistent with the predictions expected for a sample 
of normally clustered galaxies. 

(iii) The projected measurements are marginally lower 
than the predictions expected for a sample of normally clus- 
tered galaxies, and lie between the autocorrelation functions 
of normal galaxies and HII galaxies. 

Telles and Terlevich (1995) found that the space density 
of bright galaxies within IMpc 3 of HII galaxies is a factor 
~ 4 times higher than expected for a random distribution, 
but ~ 2 times less than for a control sample of sample of Sc 
galaxies. These results showed that HII galaxies are more 
clustered than a random distribution, but slightly less clus- 
tered than normal galaxies. The present work extends this 
analysis to much fainter apparent magnitudes, measuring 
the density of galaxies to Mt 7 ~ —14.5, which correspond 
to very low mass galaxies and HI clouds (~ 1O 8 M0). Our 
results are in agreement with the earlier studies, with the 
additional conclusion that HII galaxies do not have prefer- 
entially faint, low-mass neighbours. 

We conclude that star-formation in these galaxies is not 
triggered by tidal interactions and must have a different ori- 
gin, possibly associated with the formation and evolution 
of Super Stellar Clusters in starbursts, as observed in the 
UV (Vacca 1994; Meurer et al. 1995; Ho 1997) and in the 
Near-IR (Telles et al. 1999). 
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